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(57) Abstract 

The movement of a catheter (34) through a medium, 
which may be living tissue such as a human brain, is controlled 
by mechanically pushing a flexible catheter (34) having a 
magnetic tip (30) through the medium, and applying a magnetic 
field having a magnitude and a direction that guides the 
mechanically pushed catheter tip (30) stepwise along a desired 
path. The magnetic field is controlled in a magnetic stereotaxis 
system (10) by a processor (90) using an adaptation of a PID 
(proportional, integral, and derivative) feedback method. The 
magnetic fields are applied by super-conducting coils (Fig. 2), 
and the currents applied through the coils (Fig. 2) are selected 
to minimize a current metric. 
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METHOD AND APPARATUS FOR MAGNETICALLY CONTROLLING MOTION 
DIRECTION OF A MECHANICALLY PUSHED CATHETER 

BACKGROUND OF THE INVENTION 
5 1. Field of the Invention 

This invention relates to devices and methods for moving 
an implant in a body, and more particularly to such devices 
and methods that apply pushing forces with a flexible 
attachment and that magnetically steer the implant in the 
10 body with high accuracy. 

2. Description of Related Art 

There is a large body of conventional (nonmagnetic) 
stereotactic prior art, in which a frame (e.g., a so-called 

15 "BRW Frame") is attached to the skull to provide a navigation 
framework. Such a frame has arcs to determine an angle of an 
"insertion guide" which is usually a straight tube through 
which some medical therapeutic agent is passed, such as a 
biopsy tool. These methods have been confined to 

20 straightline approaches to a target. 

There is also a smaller body of prior art in which a 
handheld permanent magnet or an electromagnet is used to move 
a metallic implant. 



SUBSTITUTE SHEET (RULE 26) 



WO 99/11189 



PCT/US98/03276 



Previous implants for delivering medication or therapy 
to body tissues, and particularly brain tissue, have 
generally relied upon the navigation of tethered implants 
within vessels, or navigation of tethered or untethered 
5 implants moved intraparenchymally (in general brain tissue) 
by magnetic force. 

Navigation of untethered implants, in the past, has 
generally comprised finding ways to apply magnetic force 
optimally, including both magnitude and direction, for a 

10 given step of "free" motion. However, difficulty in finding 
a set of currents to accomplish a move step is encountered 
because of the complexity of calculating the magnetic forces 
resulting from multiple coils. 

It is well-known that two like coils on a common axis, 

15 having like currents, provide a highly uniform magnetic field 
on their axis at the midpoint between them. In addition, it 
is known that the field is approximately uniform for an 
appreciable region around the midpoint, and relatively 
strong, as compared with any other two-coil arrangement 

20 having the same coil currents. This arrangement of coils and 
currents meets the need for an accurate, strong guiding 
torgue applied to a magnetic implant near the midpoint 
between the coils. Because the field is quite uniform near 
the midpoint, undesired magnetic forces on the implant are 

25 negligible. However, this arrangement is less suitable for a 
moving implant when the implant is some distance from the 
midpoint between the coils or not on the axis, or when the 
implant axis is not along the coil axis. In these important 
cases, this simple coil arrangement cannot provide accurate 

30 directional guidance. Furthermore, simple vector 

combinations of three such coil pair arrangements cannot 
provide accurate guidance in an arbitrary direction, except 
at one spot at the center of the arrangement. 
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The Magnetic Stereotaxis System (MSS) originated from 
the hopes that a less-invasive methodology could be developed 
which would allow neurosurgeons to operate in previously 
inaccessible regions of the brain. By introducing a small 
5 permanent magnetic implant into the brain through a small 

"burr hole" drilled through the skull prior to the operation, 
large superconducting coils could be used in conjunction with 
a pushing mechanism to magnetically guide the implant and 
overlaying catheter through the brain's parenchyma, all the 

10 while avoiding the important structures of the brain. The 
operational methodology of the MSS was, and continues to be, 
expected to be less destructive to the tissues of the brain 
than the shunts, straight tubes, and other devices associated 
with conventional techniques in neurosurgery. 

15 The first MSS was conceptually developed in 1984 as the 

Video Tumor Fighter (VTF) , and is shown in U.S. Pat. No. 
4,869,247 issued Sept. 26, 1989. This system specifically 
focused on the eradication of deep-seated brain tumors via 
hyperthermia-based treatment. It was envisioned that the 

20 magnetic coils of the VTF would guide a small (-3 mm 

diameter) magnetic thermosphere through the brain into a 
tumor. Rastering the implant throughout the volume of the 
growth, the tumor cells could be destroyed by inductively 
heating the implant with radio-frequency radiation. 

25 Further studies revealed that the reality of a 

magnetomotive based system used to direct a small implant 
promised numerous applications other than the hyperthermia- 
based treatment of brain tumors by induction. These 
included: biopsy, pallidotomy, delivery of precision 

30 radiation therapy, magnetically placed implants that deliver 
chemotherapy to otherwise inaccessible tumor locations, and 
(by attaching a semi -permeable catheter to the implant) the 
delivery of chemicals to specific sites in the brain without 
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the need for penetrating the blood-brain barrier which has 
complicated contemporary systemic chemical delivery. This 
means of chemical delivery seemed particularly hopeful in the 
treatment of Parkinson's disease, where the catheter could be 
5 used to deliver dopamine to the affected regions of the brain 
with minimal indiscriminate distribution of the 
neurotransmitter to the surrounding tissue, thereby lessening 
attendant side effects. It was in the light of these 
possible broadened applications of the VTF that the system 

10 became known as the MSS . 

Referring now to Fig. 1A and Fig. IB, the most recent 
MSS apparatus 10 included six superconducting coils (not 
visible in Fig. 1A and Fig. IB) located in a rectangular box 
or helmet 12. With the z-axis defined in the direction of 

15 the axial component of the head, the x- and y- coil axes are 
rotated 45° from the sagittal plane 14 of the head, which 
would be positioned in opening 16. The x- and y- coil axes 
are symmetrically located such that the horizontal extensions 
22 of the MSS apparatus 10 away from the patient's body is 

20 minimized. Because the lower edge of the treatable part of 
the brain is typically located 10 cm above the shoulder line 
for an average adult, the z-coils {located on the body-axis 
of the supine patient) were compressed to allow for a maximum 
extension of the head into helmet 12. 

25 The vision component of the MSS consists of a 

superposition of pre-operative MRI images referenced by 
biplanar fluoroscopy cameras 20 linked to a real-time host 
system (not shown in Fig. 1A and Fig. IB) . Both cameras 20 
are calibrated to the MSS six-coil helmet design. X-ray 

30 generators for cameras 20 are located inside magnetic shields 
22. Using x-ray visible fiducial markers located on the 
skull of the conscious patient, the coordination of the 
implant's position inside the cranial volume to the helmet's 
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reference system (and hence the corresponding preoperative 
MRI scan) is done through a series of coordinate 
transformations executed by a host system and displayed for 
the surgeon on a workstation. 

The central problem to the inductively-based guidance of 
a magnetic implant pertains to the inverse problem of 
electromagnetism as influenced by Earnshaw' s theorem. The 
conventional problem of electromagnetism centers on the 
evaluation of the gradient and magnetic field given 
established magnetomotive sources. For the MSS, however, the 
situation is reversed in that the magnetic field and its 
gradient are specified at a point in space while the 
strengths of the six actuators are to be determined. Control 
of the motion and position of an untethered implant would be 
difficult in the MSS , given the fundamental instability of a 
non-diamagnetic moment in a static or quasi-static magnetic 
field as related to Earnshaw' s theorem for static/quasi- 
static magnetic fields, if it were not for the resistive 
nature of the parenchyma. In early tests, small cylindrical 
(up to 5 mm in length and 5 mm in diameter) permanently 
magnetized NdBFe objects were used. The relatively strong 
moment of these objects (0.016 A-m 2 to more than 0.04 A-m 2 ) 
facilitated the creation of the necessary aligning torque 
without the requirement of a strong magnetizing field, 
resulting in lower current values. 

The permanent magnetization of the implant requires a 
predetermined magnetic field in order to ensure that the 
implant is oriented in the desired direction. While it is 
possible to generate a magnetic force to displace the 
implant, it was found that the requirement of specific force 
and field alignment could result in unobtainable currents (as 
high as thousands of amperes) . It was also found that even 
for viable solutions, the equilibrium state was sometimes 
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unstable to such an extent that the implant tended to be 
difficult to control. 

SUMMARY OF THE INVENTION 
5 The invention is an apparatus and method for moving an 

implant in the body by applying pushing forces with a 
flexible attachment and magnetically steering the implant in 
the body with high accuracy and controllability. Because the 
intended moving force is applied non-magnetically, it is 

10 possible and desirable to apply currents in the magnetic 

steering apparatus in such combinations as to maximize the 
magnetic field at a body location inside the coil array to 
thereby provide optimal directional guidance torque on an 
implant while minimizing undesired translational force on the 

15 implant. 

According to one aspect of the invention, there is 
provided a method for controlling movement of a catheter 
through a medium, in which a flexible catheter having a 
magnetic tip is pushed through a medium, and a magnetic field 
20 having a magnitude and orientation effective to guide the 

mechanically-pushed catheter tip in a predetermined direction 
is applied. 

According to another aspect of the invention, a method 
for providing stepwise movement of a catheter having a 

25 magnetic tip is provided, in which the method includes the 
steps of selecting a desired path of the catheter through 
living tissue, inserting the catheter tip into the living 
tissue, determining actual positions of the magnetic tip and 
correction vectors (the correction vectors representing 

30 differences between locations on the desired path and the 
actual positions of the magnetic tip) , storing values of 
correction vectors in a memory, and applying a magnetic field 
adjusted to achieve movement of the magnetic tip at least 
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approximately along the desired path, the adjustment 
depending upon at least one stored set of values of 
correction vectors. 

Also provided is a device for guiding a catheter having 
5 a magnetic tip through a medium, the device comprising a 
helmet having a cavity configured to encompass a medium 
through which a catheter is to be guided, a magnetic field 
generator generating a magnetic field within the cavity, a 
position sensor sensing a location of a magnetic tip of a 

10 catheter in the cavity and generating a signal indicative of 
the sensed location, an advancement mechanism pushing the 
magnetic tip of the catheter through the medium, and a 
processor responsive to the signal from the position sensor 
and having an operator control input, the processor being 

15 configured to control the magnetic field generated by the 
magnetic field generator in response to commands input via 
the operator control input and the signal received from the 
position sensor. 

The above embodiments may also incorporate significant 

20 additional improvements, including, for example, the 

minimization of a current metric, so that the proper magnetic 
field to guide the magnetic tip through the medium is 
generated with a near-minimum amount of current. 

The methods and apparatuses of this invention provide 

25 the ability to more accurately direct a seed or catheter in 
the brain or other parts of the body, including the path to 
that position. Highly accurate directional guidance of 
implants is possible over arbitrary nonlinear paths, and the 
implant can be guided freely through tissues such as brain 

30 tissue, without being limited to the interior of vessels. 

Additional advantages of the present invention over 
prior art systems are that: 
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(1) Solutions applicable to guiding an implant on a 
predetermined path are simpler, and thus, are found more 
rapidly and with less likelihood of error for a given step of 
motion. 

5 (2) Solutions are much more stable than with prior art 

systems, and are free of runaway conditions. 

(3) Currents applied by the new method are generally 
considerably smaller than with previous methods; therefore, 
the current changes between steps are smaller, allowing 

10 changes to be made much more rapidly and accurately between 

steps, and with less possibility of guenching superconducting 
magnets . 

(4) Guidance force occurs without skid, which is a 
condition in which the magnetic field that orients the 

15 implant and the magnetic force are in different directions so 
that the axis of the implant skids along the path. 

(5) Currents are applied in a simple temporal fashion, 
moving directly from one set to another set between two steps 
of motion. The actual force impulse causing each step of 

20 motion is from the duration and distance of the externally 
applied non-magnetic force during that step. (Prior art 
systems ramped currents from conditions for subthreshold 
force to that of a moving force and then back down below 
threshold at the appropriate time, which is a complex dynamic 

25 sequence subject to substantial error in step length due to 
the tribological nature of the implant and tissue. 

(6) Navigation can now occur continuously rather than in 
steps . 

It is thus an object of the invention to provide a 
30 method for controlling the motion of a catheter in any 
predetermined direction. 
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It is a further object of the invention to control the 
motion of a catheter by applying a torque to the catheter to 
guide its direction with a reliable, predictable strength. 

It is yet another object of the invention to control the 
5 motion of a catheter rapidly, accurately, and reliably, even 
when the magnetic system used in conjunction with the 
catheter includes superconducting coils that are vulnerable 
to misoperation from too rapid current changes. 

10 BRIEF DESCRIPTION OF THE DRAWINGS 

Figs. 1A and IB are views along the Z-axis and into the 
cavity of one embodiment of a Magnetic Stereotaxis System 
(MSS) , respectively; 

Fig 2 is a simplified representation showing the 
15 orientation of superconducting coils within the helmet of an 
MSS; 

Fig. 3A is drawing of a portion of a catheter having a 
metal tip, and Fig. 3B is a drawing of a magnetic element and 
push wire within the catheter of Fig. 3A; 
20 Fig. 4 is a diagram showing various points and vectors 

along a desired path of a catheter tip, in which the points 
and vectors shown are significant for the calculation of the 
magnetic field needed to guide the tip in accordance with the 
invention; 

25 Fig. 5 is a block diagram of a portion of a Magnetic 

Stereotaxis System in accordance with the invention; 

Fig. 6 is a block diagram of another portion of the 
Magnetic Stereotaxis System of Fig. 5; and 

Fig. 7 is a block diagram of a configuration of a 
30 processor suitable for use in the MSS of Figs. 5 and 6. 
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DESCRIPTION OF THE PREFERRED EMBODIMENTS 
The invention is described herein in the context of 
guiding an implant inside a brain because tests have been 
performed for this case. However, those skilled in the art 
5 will recognize that the invention is applicable to other 

parts of the body and to other media through which a magnetic 
tip at the end of a flexible catheter can be pushed. 

In the present example, a system of six coils is 
provided. Referring to Fig. 2, which is a simplified 

10 representation of the coil geometry, the coils X+, X-, Y+, Y- 
, Z + , and Z- are operated in unbalanced pairs (X+, X-) , (Y+, 
Y-) , (Z+, Z-) acting approximately as three vectors X, Y, and 
Z having mutually perpendicular directions, but different 
magnitudes. The operation is achieved in a manner that is 

15 effective at positions off the axis of any or ail of the 
three pairs of coils (X+, X-) , (Y+, Y-) , (Z+, Z-) 
constituting the helmet (or other coil support arrangement) . 

The method for controlling the path in the Magnetic 
Stereotaxis System (MSS) includes calculations involving 

20 vectors that represent desired path steps and corrective 

feedback vectors in an inventive adaptation of a well-known 
PID (proportional, integral, and derivative) feedback method. 
However, the operations here are on the motion of a magnetic 
implant. Referring to Figs. 3A and 3B, the magnetic implant 

25 30 acts as a leading element (the "directing train engine," 

which is called a magnetic delivery vehicle or MDV) and which 
is only part of the overall moving unit which also includes a 
pushing stylette 32 driven by a motor (not shown) which 
supplies a force F. Stylette 32 and implant 30 are 

30 surrounded by flexible catheter 34, which has a metal tip 36 
at the end within which MDV 30 is located. The method and 
apparatus of this invention controls the curving elements of 
the catheter. 
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For the purposes of this discussion of the invention, 
the term "lag" is applied to the angle by which the catheter 
would fall outside a curve that would be followed by a 
completely limp catheter. During magnetic navigation, the 
MSS system monitors the position of the catheter tip (the 
MDV) versus a desired position in accordance with a planned 
path. Referring to Fig. 4, corrections are made in 
accordance with a feedback/lag model that is given by the 
equation 



in which: 

V lag is a lag-correction vector pointing to an inside of 

a local curvature of the planned path indicated by R at point 

p 2 ; and 

g is a gain factor; and 

V path is constructed by finding a point p, on the planned 
path closest to the present actual position, and then finding 
a point p 2 further along the planned path which matches the 
length of the intended step. Then V path = p, -p, ; 

V is a vector from the actual position to the point 



on the path closest to the present actual position (i.e., a 
direction opposite an "error vector"); and 

V is the resultant direction that the next step 

step 

should take, corrected for lag. 

The direction of V step is constrained to vary from move 

to move in a manner so that sudden changes in V stcp result in 
a path having a specified minimum radius of curvature or a 
less curved trajectory. The value of g is computed and 
adjusted for each step by means of a PID model in which 



step 




correction 



correction 



(l) 



correction 



SUBSTITUTE SHEET (RULE 26) 



WO 99/11189 



PCTYUS98/03276 



12 

H*,V c(n) +*,-(V c(n) + V^, +V c(B _ 2) ) + ^(V c(n) -V c(n _ 1) )j| , where : 

k p , k i , and k d are predetermined constants that may be 
found by experiment; and 

v c<n)' v c<n-i)' and V cc-2) are the correction vectors at the 
5 nth step (present step), the previous (n-l)th step, and the 
(n-2)th steps, respectively. 

It is noted that what here is called a "correction 
vector" is the negative of what is often called an "error 
vector." The term "correction vector" rather than "error 
10 vector" has been used throughout this description for 
convenience. However, this usage is not intended to be 
limiting in any way, because those skilled in the art will 
recognize the equivalence of "correction vectors" and "error 
vectors" and the notational changes that would be required by 
15 the substitution of one for the other. 

When the correction vectors at the nth, (n-l)th, and (n- 
2) th steps are used in the correction vector V step , they 
become an approximation to an integral mode of control. That 
is, a correction is made proportional to the integral of the 
20 recent path errors. Similarly, an approximate derivative or 
rate correction is made when the difference (V c(n) -V c(D _ n ) of 
the present and most recent correction vectors is used as a 
term in V step . 

In one exemplary MSS implementation, vector V lag was 
25 determined experimentally by attempting to drive a magnetic 
catheter tip (such as tip 36 in Fig. 3A) in a circle of 
predetermined radius, but setting V tag = 0 . The catheter was 
driven by a motor providing a constant advancement rate 
through a phantom gel, which provided some friction. Where 
30 there is such friction, an equilibrium arises because the 
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force from the gel counteracts the force from the magnetic 
field and tends to balance it at some point before alignment 
of movement with the magnetic field is achieved. In this 
case, the counterbalancing force of the gel simulated the 
5 counterbalancing force of actual brain tissue, and as 

expected, the magnetic tip spiraled outward from the circle 
in a direction a few degrees out of alignment with the 
magnetic field. V laf . was then experimentally set to be the 
vector correction that would have been necessary to provide 

10 to counteract the deviation from the circular path. 

By increasing the magnitude of the field, zhe magnitude 
of V, can be reduced. However, it is also desirable to at 
least approximately minimize the coil currents required to 
generate the magnetic field. Experiments have shown that, 

15 for the particular set of coils in one MSS apparatus, with a 
motor pushing a catheter along the magnetic field at a rate 
of about 0.5 mm/s, there was little decrease in V lag for 
fields greater than about 0.3 Tesla. Therefore, in the 
experiment to determine V lag , as well as in other tests, the 

20 magnetic field was constrained to have a magnitude no greate 
than 0.3 Tesla. This value has been found, by experiment, t 
be sufficient to orient the direction of movement of a 
magnetic tip of a catheter being pushed by a motor at 0.5 
mm/s through a phantom gel. It will be understood that the 

25 necessary magnetic field magnitude may vary if different 
catheters, motors, or magnetic tips are used, or if a 
different advancement rate is applied. 

Correction according to V lag is a single-parameter means 

of anticipating error that would occur due to restoring 
30 torques from an attached pushing element. A person skilled 
in the art would recognize the close relationship of this 
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action to the concept of "feed-forward" in the field of 
control theory. It is intended, according to one aspect of 
the invention, that the inventive apparatuses and techniques 
may be useful for applications in which a simple feed-forward 
5 correction is not adequate. In such applications, a computer 
or processor will have a representation of a planned or 
predetermined path stored in its memory. At any location or 
locations where the planned path deviates from a straight 
line in front of a then-present location of the seed, a 

10 program stored in the computer provides an added correction 
vector, which may be a function of several different 
parameters, to provide a better correction than V lag by 
anticipating future error in the absence of such correction. 
As one example, if a planned path curves in more than one 

15 plane, a correction vector can contain terms which have 
correction vector components in each of the planes. The 
weighting of these components preferably varies as the 
inverses of the individual radii of curvature, and are also 
preferably weighted inversely according to the distance from 

20 the present seed position at which the particular curvatures 
occur. These parameters can be readily calculated by a 
program stored in the memory of the computer. It is thus 
possible, and may be preferable, to use information related 
to a future location of the seed to guide the seed, and/or a 

25 determined rate of change in the observed correction vectors 
in addition to information stored in memory about its past 
position and errors. 

With respect to the gain parameters k p , k, , and k d , it 
has been experimentally determined that, for the tested MSS 

30 apparatus, it is sufficient to set k p =0.5, k ; =0.5, and k d =0. 
The parameter k d , which effectively adjusts the speed of 
response of the system, can be set to zero because the system 
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response obtained has been satisfactory for the relatively 
slow tip advancement rate of 0.5 mm/s. In addition, setting 
k d -0 has the advantage of reducing noise in the system, 
because the imaging method used to locate the magnetic tip 
5 (biplanar x-ray imaging) has provided an accuracy only within 
about ± 1.0 mm. Adjusting the system gain by increasing the 
magnitude of k d would result in significant noise 
amplification without a concomitant increase in position 
sensing accuracy. If more accurate methods of determining 

10 the location of the seed were used, it would be practical to 
use a nonzero value of k d to provide more rapid system 
response. This more rapid system response could permit an 
increase of the tip advancement rate above 0.5 mm/s in 
appropriate medical procedures, which may shorten the time 

15 necessary to perform such procedures. 

Once the vector V s(tp is determined, representing a 
vector corresponding to a motion of the magnetic tip at least 
approximately in the direction of the desired path, the coil 
currents necessary to generate the required field must be 

20 determined. In general, this problem is underdetermined in 
that many different sets of coil currents could be used to 
generate the required field, and, without the addition of 
further constraints, solutions can be obtained that require 
impractical amounts of currents in some or all of the coils. 

25 A method of determining practical values of coil currents 
that generate a sufficient magnitude magnetic field in the 
correct direction at the location of the magnetic tip has 
been found, and is described herein. 

For controlling superconducting or normally conducting 

30 coils, the root-mean-square value of the coil currents (the 
coil-current metric) is minimized while the strength and 
direction of the magnetic field remain at the desired values. 
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(The magnetic field is linear with respect to the coil 
currents.) This constraint removes three of the six degrees 
of freedom in selecting the six coil currents, leaving a 
quadratic equation for the current metric. This metric is 
minimized to compute the optimal coil currents. While 
minimizing the current metric does not necessarily correspond 
to minimizing the individual current magnitudes, it is a 
useful and efficient way of ensuring that the current 
magnitudes are kept at a sufficiently low level for purposes 
of the invention. 

The m equality, linearly constrained quadratic problem 
is stated as 



Maximize z = ^ xp. + X x,P u x } 

n 

Subject to X Dy Xj = e } (i=l, . . . ,m) 



;=i 



or 



Maximize z = x c + x Px( 
Subject to Dx = e 



(2) 



where xe9T. We assume that the conditions Dx = e comprise a 
non-degenerate set with m<n . If m>n , then the system is 
over specified and no solution may occur. In the case that 
m = n, then one solution exists at the point x 0 = D -1 e 
providing |D| * 0 . Hence, no optimization is possible or 
necessary . 

Constructing the Lagrangian, we find that 

L = x T c + x T Px + X T (Dx-e) (3) 

where m Lagrange multipliers have been introduced. From the 
Lagrangian, we obtain the global extremum 

(4) 

where it is assumed that the matrix inversion is possible. 

To particularize this result to the minimization of the 
current metric for the MSS, we begin by focusing solely on 





2P D T " 






C:H 




■ - 






D 0 
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the static and quasi-static cases, thus assuming that the 
source currents are either held constant or ramped relatively 
slowly in time. Given n magnetomotive sources (or r in this 
sense, actuators), we wish to operate the sources with 
5 minimal currents so that a desired magnetic field may be 

specified for a selected point in space x 0 where x„ e<R . The 
total magnetic field at any point x r b(x) , is the linear 
superposition of the magnetic fields due to each source 
evaluated at x: 

10 b(x) = Xb,(x) (5) 

M 

Since b,.(x) is linear with respect to its corresponding source 
current,/,, the above may be written as 

b(x,I) = XB f (x)/,. (6) 

where B,(x) consists of the three current-independent 
15 components of the magnetic field for each source. If we 
define the 3x« matrix B(x) as 

B(x) = {B,(x) I B 2 (x) | ••• | B„(x)} (7) 
and we write the currents as the n-element column vector I, 
then Eq. (6) and Eq. (7) can be combined to form the matrix 
20 relationship 

b(x,I) = B(x)I (8) 
Note that «>3 is assumed in order for an optimization to be 
made for a desired magnetic field b (which consists of three 
components). If «<3 (i.e., two actuators or less), the 
25 system is over constrained and no solution exists unless 

there is a degeneracy in Eq. (8). If n = 3, then the solution 

to the currents, I 0 , is given by I 0 = B"'(x 0 )b . 

We now focus our attention on the current metric defined 

as 
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(10) 



z(l) = £i? = I T I (9) 
i=i 

While the metric does not specifically limit the individual 
currents, it does serve as a means of penalizing those 
solutions that are associated with strong currents. The 
5 problem of finding an optimal set of currents (for n>3 
sources) may now be stated in a form for which there are 
7M = 3 equality constraints: 
Maximize z(l) = -I T I 

> 

Subject to B(x 0 )l = b 

It is noted that we could just as easily minimize the current 
10 metric above without loss of generality; however, writing the 
problem in the form of Eq. (10) is convenient for our present 
purposes. Since the metric of Eq. (10) is a concave 
function, the solution for the currents is: 

'-2I„ B(x.)TjO] 

>W o J I 

15 where I D is the identity matrix. Inspection of Eq. (10) 
reveals that only one extremum (and hence, the global 
maximum) of the negative current metric can occur. Eq. (11) 
thus provides the solution for the coil currents I 0 in the MSS 

given a specified magnetic field. 

20 While the current metric is sufficient in restricting 

the currents to small values in most cases, it does not 
minimize them. It is possible that a larger metric results 
when smaller currents are distributed over several sources. 
For example, say the desired field for a four source system 

25 corresponds to the optimal set of currents 

I o T = {l0 80 10 80}(A) for which z(l 0 ) = 13000 A 2 . If the 
individual currents must be less than 75 A, another (possibly 
more useful solution in some cases) would correspond to 



-i 

(0| 

(11) 
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I o T = {60 70 60 70}(A), for which z(l 0 ) = 17000 A 2 , providing the 
currents generate the same magnetic field. Including the k 
linear current limits DI>e into Eq. (10), our general n- 
source, linearly constrained problem is stated as 
Maximize z = -I I 

Subject to B(x 0 )l = b (12) 
DI>e 

Since it is more commonly found that the n actuators possess 
upper and lower limits according to |/,.|</ m „ {i -\,...,n) , the 
constraints form a closed and bounded set providing the 
specification of B still holds for the range of allowed 
currents. The problem becomes 



Maximize z(l)=-I T I 
Subject to B~(x 0 )l = b 




(13) 



where k = 2n inequality constraints with I maxJ = I nax for i = \,...,n 
have been introduced. The conditions that must be satisfied 
in order for maximum to exist are given by 

-2i 0+ u: M+*;B(x o ) = 0 



[ 


'-Id 1 




^max 


1) 








I™, 





-0 



B(x 0 )I 0 =b 
-I ± 

•■■max 



-I <I <I 

A max — *o — max 



(14) 



The possible 2 2 "=4" solutions of the above set of equations 
follow from Eq. 2 9 of Appendix B where individual constraints 
are activated among the In inequality conditions. As was 



SUBSTITUTE SHEET (RULE 26) 



WO 99/1 1 1 89 PCT/US98/03276 

20 

previously discussed, when the activated constraints combined 
with the equality constraints outnumber the degrees of 
freedom, the system of equations become over specified and no 
solution need be calculated. For those cases in which the 
5 system of equations is exactly specified, the solution must 
be checked against the inequality constraints to deem it 

viable. There remain 4" - V 7 \ — of the 4" cases which 

£ 3 (2«-/)M 

can be solved (assuming a solution exists) . Those solutions 
that satisfy the constraints are saved and the set that 

10 results in the maximum value of z(l o ) = -I 0 T I o [or minimum of 
z(l 0 ) = I 0 T I 0 ] is reported as the optimal solution. 



It is sometimes useful to restrict the magnetically 
generated force on a small permanent moment. For example, 
15 quasi-static systems such as magnetic suspensions and the MSS 
can profit from an inclusion of force constraints if higher 

currents are acceptable. The force at x o , f(x o ) , generated on 
a small permanent magnetic moment, m, due to a magnetic 
field, b, is given by 

20 f(x 0 ) = v(m T b(x))|^ (15) 

An easier notation for the present purposes involves writing 
the three component of the force as 

/,(«.)— (. = 1,2,3) (16) 

For those problems in which the moment is allowed to align 

|| m | 

25 with the magnetic field [i.e., m = "jy b ( x o) ] / Eq. (16) is 



transformed into 



\r( db(x) 

dX; j 



(z = l,2,3) (17) 
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where the strength of the moment ||m|| is known. Combining the 
results of Eq. (17) with Eq. (8), a somewhat complicated 
problem arises for those cases in which the orientation of 
the magnetic field at x 0 is unrestricted. This can be seen 
in the nonlinear form of 

For the present purposes, only those cases that rely on a 
specified magnetic field are considered. 

Using the current dependence of the magneric sources and 

a predetermined magnetic field b where b = B(x o )l, Eq. (17) can 
be written in two forms. The linear and quadratic forms are 
given by, respectively, 

'SB 



and 



While the form of Eq. (19) may appear more useful, at least 
seven actuators must be present in order to overcome the six 
constraints due to the specification of the magnetic field 
and force. If it is important that both the force and field 

20 be specified and if there are a sufficient number of 

actuators, then the work follows from Eq. (10) - Eq. (14) 
with the three additional force constraints being included 
into the field constraints. If there are exactly six 
actuators, then there exists a unique solution to the problem 

25 at 
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I 0=1 l| b f V ^ (x) t 
I B 



(21) 



providing the operating matrix is invertible and the currents 
are unbounded. 

More often than not, the experimenter is more concerned 
with either minimizing a component of the force or the 
strength of the force with respect to a limited range of 
currents and a desired magnetic field rather than specifying 
a specific value of the force. If examining a component of 
the force, Eq. (20) is generalized so that the force along 
the unit vector u (||u|| = l) is mxnimized. The force component 

of interest becomes u T f(x 0 ,l) and the problem is written as 
Maximize z(l) = -u T f(x 0 ,l) = £ u, I T B(x 0 ) T f||-(x 0 )l I 



Subject to B(x 0 )l = b 



\;=1.2.3 



-Ir 



I>- 



(22) 



Likewise, if the force strength is to be minimized, a 
quadratic form is obtained for the objective function by 
squaring the force components of Eq. (19) : 



Maxbmze z(l) = -f(x 0> l) T f(x 0 ,l) = -[=J [f^J>b^(x 0 )J I 

Subject to B~(x 0 )l = b 



SB 



f 1.1 






H 











(23) 



If it is desired that the force be maximized rather than 
minimized, then the negative sign to the objective function 
is left off in Eq. (22) and Eq. (23) . The conditions that 
establish the existence of a minimum or maximum follow from 
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Eq. 27 of Appendix A. Note that only for force minimizations 
may the currents be left unbounded. 

Portions of the above description of hardware and 
control methods and apparatuses refer to a six-coil system 
5 for use in certain applications in the medical field. 

However, the invention is not limited to systems having any 
particular number of coils, and is suitable for applications 
(including magnetic surgery applications) for which other 
coil numbers and arrangements may be preferable. For 

10 example, it is possible to remove the front coil of the above 
system, and to bring side coils closer together. The 
resulting five-coil arrangement would allow the patient to be 
more accessible to a surgeon, an important consideration for 
some surgical procedures. Other arrangements with different 

15 numbers of coils may be particularly useful with some other 
types of operations, for example, those at locations in the 
body other than the head. (It should be noted that it is not 
required by the invention that embodiments of multiple coil 
systems necessarily have coils that are operated in pairs, 

20 opposed or otherwise. Thus, the invention also encompasses 
systems having arbitrary numbers of coils.) 

An apparatus that controls coil currents and magnetic 
tip advancement is illustrated in block diagram form in Fig. 
5 and Fig. 6. (It should be noted that not all external 

25 connections to console 54 are shown in each figure.) Fig. 5 
shows a catheter 34 with metal tip 36 having a magnetic 
implant 30 inside (see also Fig. 3A) . This catheter is shown 
with tip 36 already implanted in the brain of a patient 50, 
at the beginning stages of a medical procedure. Push wire 32 

30 is operatively coupled to a motor in an advancement mechanism 
52 that a surgeon may control by issuing commands at a 
console 54. Console 54 includes a real-time computer having 
a display screen 56 and an operator control 58, which may be 
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a keypad or other convenient data entry device. Although the 
apparatus is described in conjunction with operations on a 
catheter in the brain of patient 50, it should be recognized 
that the inventive apparatus and techniques may be applied to 
5 other living tissues as well, or in other media, living or 
not, through which it may be desired to push a magnetic tip 
on a flexible push wire or guide wire. 

Referring now to Fig. 6., the head of patient 50 (not 
shown in Fig. 6) is disposed in opening 16 of coil apparatus 

10 10, as the surgeon operates the inventive apparatus. Console 
54 contains a processor such as a digital computer responsive 
to operator commands input through operator control 58 and to 
cameras 20 for controlling a power supply (not shown) that 
supplies currents to coil terminals 24 to generate the 

15 required magnetic fields. To fully automate the catheter 

movement, advancement mechanism 52 may also be controlled by 
the processor. Cameras 20 and X-ray generators inside 
magnetic shields 22 provide magnetic tip position information 
to the processor in console 54. It will be recognized, 

20 however, that other suitable means for providing position 

information may be substituted for the cameras 20 and X-ray 
generators . 

To perform a procedure, a surgeon would observe the 
initial position of the magnetic tip 36 with the aid of 

25 console 54 and plan a path through the tissue of patient 50 

(in this case, the brain) to a point at which treatment is to 
be delivered. The surgeon may also choose a step size and an 
advancement rate, although either or both of these may be 
preprogrammed. The surgeon then activates the system by 

30 issuing appropriate commands through console 54. 

In response to these commands, the processor inside 
console 54 computes positions of metal tip 36 from data 
received from cameras 20, as well as vectors V and the 
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coil currents necessary to move metal tip 36 along the 
desired path at the beginning of each step, and applies the 
computed currents to the coils while advancement mechanism 52 
advances metal tip 36 through the tissue. Correction vectors 
5 V are stored in memory to provide the necessary history 

correction J *■ 

for the PID calculations. The advancement continues until 
the endpoint of the path is reached. Although the 
advancement occurs in steps, the positions and coil currents 
can typically be calculated rapidly enough to make the 

10 advancement of magnetic tip 36 essentially continuous. 

If necessary, the surgeon may intervene by stopping 
movement of the catheter or by changing its course before the 
originally planned endpoint is reached. When a course change 
is made, the correction vectors V correctjon stored in the memory 

15 of the processor during advancement on the aborted path are 
either cleared from memory or simply ignored so that these 
vectors do not influence the calculations performed to direct 
the magnetic tip along the new path. 

Fig. 7 shows a simple block diagram of a manner in which 

20 a processor 90 may be configured for use in console 54. 

Processor 90 is provided with operator control 58 (which may, 
for example, be a keyboard) for accepting operator input and 
a display 56 for displaying information to the operator. 
Inputs and outputs 70, 74, 76, 78, and 80 corresponding to 

25 those shown in Fig. 5 and Fig. 6 are provided to obtain data 
required for sensing position of the catheter tip and for 
controlling the X-ray fluoroscopes and superconducting coils 
inside helmet 12. A random access memory 94 is provided for 
storing values of V correction , and a programmable power supply 92 

30 is provided to supply current to the superconducting coils 
through a plurality of lines 72 in accordance with the 
current values computed by processor 90. Although 
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programmable power supply 92 is shown as a power supply 
having multiple, independent, programmable outputs, it will 
be readily appreciated that a plurality of individual power 
supplies may be used as an alternative to the power supply 
5 shown in the figure. 

A document describing the structure of a computer 
program to operate a processor in accordance with the 
invention is attached as Appendix A. Computer programs 
written in the C++ computer language that determine and 

10 control current applied to the coils in accordance with the 
invention appear in Appendix B. A more detailed treatment of 
guadratic optimizations pertaining to current solutions for 
magnetic field sources appears in Appendix C. A more 
detailed treatment of the generation of a magnetic field for 

15 a circular coil of distributed current appears in Appendix D. 
A detailed mathematical treatment of the summation of field 
components for single and multiple coils appears in Appendix 
E, which continues the discussion of the material in Appendix 
D. 

20 The above-described embodiments are intended to be 

illustrative only. For example, there are numerous types of 
magnetic surgery procedures for which the coil systems 
described and the method of calculating and providing 
currents to generate fields and forces are important, but for 

25 which there is no planned or predetermined path and no 

feedback. Instead, a device in accordance with the general 
principles of the invention can provide magnetic guidance of 
a tip of a surgical device such as a catheter, endoscope, 
etc. The invention can be readily adapted so that a surgeon, 

30 under guidance from an imaging system, uses the magnetic 

system to negotiate otherwise difficult turns and movements 
of the surgical device as he or she pushes a device along the 
interior of a vessel. It will also be recognized that many 
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of the inventive methods and apparatuses may be used in 
conjunction with any coil in a non-resonant circuit that 
applies a magnetic force on a suspended or embedded object 
that is magnetically moveable. Many other modifications 
5 falling within the spirit of the invention will be apparent 
to those skilled in the art. Therefore, the scope of the 
invention should be determined by reference to the claims 
below and the full range of equivalents in accordance with 
applicable law. 
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WHAT IS CLAIMED IS: 

1. A method for controlling movement of a catheter 
through a medium, the method comprising the steps of: 

(a) mechanically pushing a flexible catheter having a 
magnetic tip through the medium; and 
5 (b) applying a magnetic field having a magnitude and an 

orientation effective to guide the mechanically-pushed 
catheter tip in a predetermined direction. 

2. The method of claim 1 and further comprising the 
steps of sensing catheter tip positions, and adjusting the 
applied magnetic field so that the catheter tip advances 
stepwise at least approximately along a predetermined path. 

5 

3. The method of claim 1 and further comprising the 

steps of: 

sensing catheter tip positions; 

generating a signal indicative of the position of the 
5 catheter tip; 

determining a correction vector between the position of 
the catheter tip and the predetermined path at least 
partially in response to said generated signal; and 

adjusting the applied magnetic field at least partially 
10 in response to the determined correction vector so that the 
catheter tip advances at least approximately along a 
predetermined path. 

4. The method of claim 3 further comprising the step of 
storing correction vectors in a memory, and wherein the 
determining step comprises determining the adjustment at 
least partially in response to a presently-determined 
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5 correction vector and at least one of the previously stored 
correction vectors. 

5. The method of claim 1 wherein the step of applying a 
magnetic field comprises the step of applying currents to a 
plurality of superconducting coils. 

6. The method of claim 5 wherein the currents are 
selected to minimize a predetermined current metric. 

7. The method of claim 5 wherein the currents are 
selected to minimize a current-limiting constraint set. 

8. The method of claim 1 wherein the step of 
mechanically pushing the catheter comprises the step of 
pushing the catheter through living tissue. 

9. The method of claim 8 wherein the living tissue is a 
human brain. 

10. The method of claim 8 further comprising the steps 
of sensing catheter tip positions, and adjusting the applied 
magnetic field so that the catheter tip advances stepwise at 
least approximately along a predetermined path. 

5 

11. The method of claim 10 further comprising the steps 

of: 

generating a signal indicative of the position of the 
catheter tip; 

5 determining a correction vector between the position of 

the catheter tip and the predetermined path at least 
partially in response to said signal; 
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and wherein the adjustment step comprises determining 
the adjustment to the applied magnetic field at least 
10 partially in response to the determined correction vector. 

12. The method of claim 11 further comprising the step 
of storing correction vectors in a memory, and wherein the 
determining step comprises determining the adjustment at 
least partially in response to a presently-determined 

5 correction vector and at least one previously stored 
correction vector. 

13. The method of claim 11 further comprising the step 
of storing correction vectors in a memory, and wherein the 
determining step comprises determining the adjustment at 
least partially in response to a rate of change of correction 

5 vectors in a previous step. 

14. The method of claim 11 further comprising the step 
of storing correction vectors in a memory, and wherein the 
determining step comprises determining the adjustment at 
least partially in response to an anticipated future error 

5 with respect to a planned path. 

15. The method of claim 8 wherein the step of applying 
a magnetic field comprises the step of applying currents to a 
plurality of superconducting coils. 

16. The method of claim 15 wherein the currents are 
selected to minimize a predetermined current metric. 

17. The method of claim 15 wherein the currents are 
selected to minimize a current-limiting constraint set. 
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18. A method for stepwise movement of a catheter having 
a magnetic tip in a magnetic field, the method comprising the 
steps of: 

(a) selecting a desired path of the catheter through 
5 living tissue; 

(b) inserting the catheter tip into the living tissue; 

(c) determining actual positions of the magnetic tip and 
associated correction vectors, the correction vectors 
representing differences between locations on the desired 

10 path and the actual positions of the magnetic tip; 

(d) storing said correction vectors in a memory; and 

(e) applying a magnetic field adjusted to achieve 
movement of the magnetic tip at least approximately along the 
desired path, the adjustment depending upon at least one 

15 stored correction vector. 

19. The method of claim 18 wherein the desired path has 
an endpoint, and the applying step is performed a plurality 
of times until the magnetic tip approximately reaches the 
endpoint . 

5 

20. The method of claim 18 wherein the determining step 
comprises determining a point pi on the desired path closest 
to the then actual position of the magnetic tip and computing 
a vector V correctjon representing the distance and correction from 

5 the then actual position of the magnetic tip to the point pi. 

21. The method of claim 18 wherein the applying step 
comprises applying the magnetic field in accordance with a 
predetermined function of at least a present correction 
vector, a previously stored correction vector, a lag- 

5 correction vector, and a vector along the desired path having 
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a predetermined length, so that the magnetic tip advances at 
least approximately along the desired path. 

22. The method of claim 21 wherein the determining step 
comprises determining a point pi on the desired path closest 
to the then actual position of the magnetic tip and computing 
a vector V correction representing the distance and correction from 

5 the then actual position of the magnetic tip to the point pi. 

23. The method of claim 22 wherein the step of applying 
the magnetic field comprises: 

selecting a point p 2 a predetermined distance from pi 
along the desired path; 
5 calculating a vector V path directed from pi to p 2 ; and 

changing the magnetic field to move the magnetic tip in 
accordance with a function of a vector V sl4p , where 

V 

V, =V ,,,+V, +g c ° rrecti<>n , in which: 

T step * path lag 5 i y i ' 

I correction I 

V lag is a lag-correction vector pointing to an inside of 

10 a local curvature of the planned path at point p 2 ; and 
j is a gain factor. 



24. The method of claim 18 wherein the applying step is 
performed a plurality of times, each time resulting in 
movement of the magnetic tip from a then current position to 
another position, and at least some of the applying steps are 

5 dependent upon the correction vectors determined for 

positions of the magnetic tip prior to the then current 
position . 

25. The method of claim 24 wherein the step of storing 
correction vectors comprises retaining at least two 
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consecutively determined correction vectors in memory, each 
determined for positions of the magnetic tip prior to the 
5 then current position, and the applied magnetic field is 

dependent upon a function of at least the two consecutively 
determined correction vectors. 



26. The method of claim 25 wherein: 

(a) the step of determining the correction vectors 
comprises, for each correction vector: 

determining a point p x on the desired path closest 
5 to the then current position of the magnetic tip; and 

determining a correction vector V correction representing 
the distance and direction from the then current position of 
the magnetic tip to the point pi; and 

(b) the step of applying the magnetic field further 
10 comprises: 

selecting a point p 2 a predetermined distance 
further along the desired path from pi; 

calculating a vector V path representing the distance 

and direction from pi to p 2 ; 
15 changing the magnetic field to move the magnetic 

tip in accordance with a function of a vector V 5tep , where 
V 

V =V «.+V. +e correct ""' , and where: 

step path la S * I VT I 

1 correction I 

V lag is a lag-correction vector pointing to an 
inside of a local curvature of the planned path at point pi; 
20 HM«<-> + W») + V «-> +V c(n- 2 )) + M V c ( n, "V c(n _ 0 )| , where 

k , k, , and k d are predetermined constants; and 

V c(n) , V B( ._ 1)f and V c(n _ 2) are a correction vector 
determined for the then current position of the magnetic tip, 
a correction vector determined for a position of the magnetic 
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tip before an immediately preceding magnetic field 
application, and a correction vector determined for a 
position of the magnetic tip before a second immediately 
preceding magnetic field application, respectively. 

27. The method of claim 18 further comprising selecting 
a minimum radius of curvature for movement of the magnetic 
tip; and wherein the step of applying the magnetic field 
comprises constraining said magnetic field so that the 
magnetic tip follows a trajectory having a curvature not less 
than a predetermined minimum radius of curvature. 

28. A method for moving a catheter having a magnetic 
tip disposed in a magnetic field stepwise along a path, the 
method comprising performing the steps of claim 18 first for 
a first desired path and then for a second desired path, and 
wherein applications of the magnetic field are responsive 
solely to correction vectors determined with respect to the 
respective paths. 

29. The method of claim 18 wherein the step of 
selecting a desired path comprises: 

selecting an ideal path; 

approximating the ideal path by parametric cubic 
polynomials, whereby the approximated path becomes the 
desired path; and 

storing a representation of the parametric cubic 
polynomials in a memory. 

30. The method of claim 29 further comprising selecting 
a minimum radius of curvature for movement of the magnetic 
tip; and wherein the step of applying the magnetic field 
comprises constraining said magnetic field so that the 
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magnetic tip follows a trajectory having a curvature not less 
than a predetermined minimum radius of curvature. 

31. The method of claim 18 wherein the application step 
comprises controlling currents through a plurality of 
superconducting coils. 

32. The method of claim 18 wherein the application step 
comprises controlling currents through a plurality of coils, 
wherein the currents in the coils are characterized by a root 
mean sguare value, and further wherein the application step 
includes minimizing the root mean square value of the coil 
currents for each application. 

33. The method of claim 32 wherein the application step 
comprises adjusting currents in three pairs of coils. 

34. A device for guiding a catheter having a magnetic 
tip through a medium, the device comprising: 

(a) a helmet having a cavity configured to encompass a 
medium through which a catheter is to be guided; 

(b) a magnetic field generator for generating a magnetic 
field in the cavity; 

(c) a position sensor for sensing a location of a 
magnetic tip of a catheter in the cavity and generating a 
signal indicative of the sensed location; 

(d) an advancement mechanism for pushing the magnetic 
tip of the catheter through the medium; and 

(e) a processor responsive to the signal from the 
position sensor and having an operator control input, the 
processor being configured to control the magnetic field 
generated by the magnetic field generator in response to 
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commands input via the operator control input and the signal 
received from the position sensor. 

35. The device of claim 34 wherein the magnetic field 
generator comprises a set of six superconducting coils. 

36. The device of claim 35 wherein the six 
superconducting coils are configured as three sets of 
unbalanced pairs generating magnetic fields in mutually 
perpendicular directions. 

37. The device of claim 36 wherein the helmet is 
configured to encompass a human head. 

38. The device of claim 34 wherein the magnetic field 
generator comprises a plurality of coils. 
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